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Abstract 

The flavor singlet axial charge has been a source of study in the last years due to its relation 
to the so called Proton Spin Problem. The relevant flavor singlet axial current is anomalous, 
i.e., its divergence contains a piece which is the celebrated Ua(X) anomaly. This anomaly is 
intimately associated with the rj meson, which gets its mass from it. When the gauge degrees of 
freedom of QCD are confined within a volume as is presently understood, the Ua{^-) anomaly is 
known to induce color anomaly leading to "leakage" of the color out of the confined volume (or 
bag). For consistency of the theory, this anomaly should be canceled by a boundary term. This 
"color boundary term" inherits part or most of the dynamics of the volume (i.e., QCD). In this 
paper, we exploit this mapping of the volume to the surafce via the color boundary condition to 
perform a complete analysis of the flavor singlet axial charge in the chiral bag model using the 
Cheshire Cat Principle. This enables us to obtain the hitherto missing piece in the axial charge 
associated with the gluon Casimir energies. The result is that the flavor singlet axial charge is 
small independent of the confinement (bag) size ranging from the skyrmion picture to the MIT 
bag picture, thereby confirming the (albeit approximate) Cheshire Cat phenomenon. . 



Pacs: 12.39-x, 13.60.Hb, 14.65-q, 14.70Dj 
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1 Introduction 



The possibility of formulating a physical theory by means of equivalent field theories defined in 
terms of different field variables, leads to a construction principle for phenomenologically sen- 
sible and conceptually powerful models, referred to as the Chesire Cat Principle (CCP)|I], ||]. 
In 1+1 dimensions fermionic theories are bosonizable Q and the CCP can be made exact and 
transparent. Furthermore in the supersymmetric world, the powerful nonrenormalizable theo- 
rems allow a web of exactly equivalent theories to be established. In the real four-dimensional 
nonsupersymmetric world, bosonization with a finite number of degrees of freedom is not exact. 
However based on the unproven "theorem" of Weinberg Q, it seems possible to argue that the 
CCP should hold also in four dimensions, albeit approximately. In view of recent developments 
on establishing the network of dualities in what is believed to be a fundamental theory (i.e., 
string theory) where one might say that the exact CCP holds, the notion of a precise CCP in 
the real world is no longer so preposterous. The aim of the present investigation is to show the 
full consistency of the CCP in the hadronic world for the case of the Proton Spin, which was not 
fully satisfactorily established in our previous efforts in this direction |||, 0Q. We complete 
the program in this paper. 

Quantum Chromodynamics (QCD) is the theory of the hadronic phenomena Q. At suffi- 
ciently low energies or long distances and for a large number of colors Nc, it can be described 
accurately by an effective field theory in terms of meson fields ||. In this regime, the color 
fermionic description of the theory is extremely complex due to confinement. However the im- 
plementation of the CCP in a two phase scenario called the Chiral Bag Model (CBM) has proven 
surprisingly powerful [|l(| . 

What is the CBM? Let space-time be divided in two regions by a hypertube, that is, the 
evolving bag. In the interior of the tube, the dynamics is defined in terms of the microscopic 
QCD degrees of freedom, quarks and gluons. In the exterior, one assumes an equivalent dynamics 
in terms of meson fields, i.e., one that respects the symmetries of the original theory and the 
basic postulates of quantum field theory The two descriptions are matched by defining 
the appropriate boundary conditions which implement the symmetries and confinement |L], [lC| ]. 
What this does effectively is to delegate all or part of the principal elements of the dynamics 
taking place inside (QCD) the bag to the boundary. We will see that this strategy works quite 
efficiently in the problem at hand. 

In this scenario the CCP states that the hadron physics should be approximately independent 
of the spatial size of the confinement region or the bag Q . This realization of the principle has 
been tested in many instances in hadronic physics with fair success ]2|. 

There is one case, however, where the realization of the CCP has not been as successful as in 
the other cases, namely, the calculation of the flavor singlet axial charge (FSAC) of the nucleon. 
Indeed in the previous efforts |5|, [7|], the CCP was realized only partially as it seemed to fail 
at certain points such as for zero bag radius. It is the leitmotiv of this work to remove this 
apparent failure. 

The observable FSAC has become very relevant in the nucleon structure in recent years, 
because it is associated with the so called Proton Spin Problem |ll]]. The experimentally observed 
small value for the FSAC implies a strong violation of the so-called Ellis-Jaffe sum rule and 

1 Note that in these papers, we have shown that the CCP holds for non-zero bag radii but it failed when the 
bag radius shrank to a point, implying that in the model studied, the pure skyrmion and the MIT bag did not 
have the equivalent structure required by the CCP. 
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therefore implies that the polarization of the proton is not carried exclusively by the valence 
quarks. It is also very interesting from the formal point of view, because the flavor singlet axial 
current - which is the origin of the observable - is anomalous, and its anomaly is related in a 



non-trivial way with the gluonic structure of the theory [13]. 

In the CBM, the scenario of how the CCP is realized - which is the central issue of our problem 
- is very intricate. As stated, the flavor singlet axial current is associated with the anomaly and 
effectively with the 7/ meson. Thus, besides the pion field of the conventional effective theories 
which accounts for spontaneously broken chiral symmetry, the correct treatment of the FSAC 
requires minimally the inclusion of a field describing the rj meson. We shall label it for brevity 
tj(x) since no confusion will arise in what follows. 

The intricacies of the hedgehog configuration and its relevance to the fractionation of baryon 
charge and other observables have been extensively discussed [Q] and fairly well understood 



15, They will be implemented in our calculation without much details. Moreover the 

inclusion of the rj meson carries subtleties of its own. The vacuum fluctuations inside the bag, 
that induce the baryon number leakage into the skyrmion |jl4|], also induce a color leakage if a 
coupling to a pseudoscalar isoscalar field is allowed [fT?f . This leakage would break color gauge 
invariance and confinement in the model unless it is canceled. As suggested in [O], this color 
leakage can be prevented by introducing into the CBM Lagrangian a counter term of the form 

C C t = j> df3K^{TAnU^ - TrlnC/) (1) 

where Np is the number of flavors (here taken to be =3), j3 is a point on a surface E, is the 
outward normal to the bag surface, U is the U(Np) matrix- valued field written as U = e l?r// ^ e lT, l f 
and the properly regularized Chern-Simons current = e^ va P (G^G^p — \f abc gG^G b a Gp) 
given in terms of the color gauge field G^. Note that (||) manifestly breaks color gauge invariance 
(both large and small, the latter due to the bag), so the action of the chiral bag model with this 
term is not gauge invariant at the classical level but as shown in [O] , when quantum fluctuations 
are calculated, there appears an induced anomaly term on the surface which exactly cancels this 
term. Thus gauge invariance is restored at the quantum level. 

The equations of motion for the gluon and quark fields inside and the rj field outside are the 
same as in [||, |6[. However the boundary conditions on the surface with the inclusion of Eq.(|l|) 
read @ 

n-E a = ^-h ■ B a n (2) 

8vr 2 / ' V ; 

Nfq 2 

hxB a = -^-h x E a r, (3) 

and 

1 - N F g 2 

-n • (^775^) = fh-dr]+ ^r^~™ • K (4) 

where E a and B a are, respectively, the color electric and color magnetic fields. Here ip is the 
QCD quark field. 

A complete treatment calls for a full Casimir calculation of the gluon modes, which is highly 
subtle due to the p-wave structure of the r\- field. Such a calculation is in progress |18| and 
will be reported in a later publication. Here we would like to side-step this technically difficult 
procedure by first assuming the CCP in evaluating the Casimir contribution with the color 
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boundary conditions (||), (|3|) and (Q) taken into account and check a posteriori that there is 
consistency between the assumption and the result. 

The next section will define our formulation, recall our old results, and clarify the new 
contributions. Section 3 will focus on the gluon Casimir contribution to the FSAC, the major 
contribution of this presentation. Finally section 4 will contain the results, conclusions and the 
future prospects of ongoing work. 



2 The Chiral Bag Formalism 

Our aim is to calculate the FSAC in the CBM scenario. In order to do so we need a spe- 
cific formulation of the model through its equations of motion and boundary conditions. The 
equations of motion have been shown repeatedly in our previous works |6L [7l and the color 
boundary conditions were recalled in the introduction. We refer the reader to those references 
for a detailed discussion on their structure, their resolution and the implementation of gauge 
invariance and confinement. Our calculation will be carried out in the static spherical cavity 
approximation, that is, our bag will be a static sphere of radius R dividing two regions of space 
in which the theory is implemented by QCD for r < R, and by an effective meson theory for 
r > R. 



2.1 The anomaly and proton spin 

The anomalous suppression of the first moment, T^, of the polarized proton structure function 
g\ has been the focus of intense theoretical and experimental activity for nearly a decade. While 
it is now generally accepted that the key to understanding this effect is the existence of the 
chiral U(l) anomaly in the flavor singlet axial current there are several explanations reflecting 
different theoretical approaches to proton structure. 

The starting point is the sum rule for the first moment, i.e., 

Tf(Q 2 ) ^ £ dxgf&Q 2 ) = ±C? s (a s (Q 2 )) (a 3 + ^a 8 ) 

+ icf(a 5 (Q 2 ))a°(Q 2 ). (5) 

Here Ci(a s ) are first moments of the Wilson coefficients of the the singlet (S) and nonsinglet 
(NS) axial currents and a s the perturbatively running QCD coupling constant. Moreover a 3 , 
a 8 and a°(Q 2 ) are the form factors in the forward proton matrix elements of the renormalized 
axial current, i.e., 

(P, s\A^\p, s) = s^a 3 , (p, slA^p, s) = Sfl ^=a 8 , 

and 

(v,s\A^\p,s) = s^a°, 

where p^ and are the momentum and the polarization vector of the proton, a 3 and a 8 can 
be chosen Q 2 independent and may be determined from the and ratios. a°(Q 2 ) evolves 
due to the anomaly and its evolution can be described in the AB scheme [11| by 

a (Q 2 ) = AZ-N F ^^Ag(Q 2 ). (6) 
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Naive models or the OZI approximation to QCD lead at low energies to 

a ps a 8 ps 0.69 ± 0.06. (7) 



Experimentally [ 19 1 



a°(oo) = 0.10i^ 7 . (8) 
The explanation for this unexpected small value has given rise to many interpretations related 



to hadron structure, vacuum structure and evolution [20, 21, 22 



2.2 The formalism 

To obtain the FSAC, we need to calculate the matrix elements of the flavor singlet axial current. 
Let us write the current in the CBM as a sum of two terms, one from the interior of the bag 
and the other from the outside populated by the meson field rj (we will ignore the Goldstone 
pion fields for the moment; they will be taken into account for the baryon charge leakage) 

A" = A B @ B + A%Q M . (9) 

Since we will be dealing only with the flavor-singlet axial current, we will omit the flavor index 
in the current. We shall use the short-hand notations 0_b = 9{R — r) and Qm = 9{ r — R) with 
R being the radius of the bag. We demand that the Ua{X) anomaly be given in this model by 

Q ^ = ^ ~ a + fm 2 v@M (1Q) 

a 

Our task is to construct the FSAC in the chiral bag model that is gauge-invariant and consistent 
with this anomaly equation. Our basic assumption is that in the nonperturbative sector outside 
of the bag, the only relevant £^(1) degree of freedom is the massive rj field. This assumption 
allows us to write 

A% = A» = fd»r, (11) 

with the divergence 

= fm 2 vV . (12) 

Now the question is: what is the gauge-invariant and regularized A B such that the anomaly ( |lO| ) 
is satisfied? To address this question, we rewrite the current (||) as 

A^ = A" Bq +A Bg +A^ (13) 

such that 

dp(A* Q +A%) = fm 2 v r,e M , (14) 

9/b g = ^££°-£°e s . (15) 

a 

The subindices Q and G imply that these currents are written in terms of quark and gluon fields 
respectively. In writing (|l4|), we have ignored the up and down quark masses. We should stress 
that since we are dealing with an interacting theory, there is no unique way to separate the 
different contributions from the gluon, quark and rj components. In particular, the separation 
we adopt, dl4j ) and (|15D, is non- unique although the sum is without ambiguity . We found 
however that this separation leads to a natural partition of the contributions in the framework 
of the bag description for the confinement mechanism that we are using here. 
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Figure 1: Various contributions to the flavor singlet axial current of the proton as a function of 
bag radius : (a) quark contribution Ab q ', (b) 7/ contribution A v and (c) the sum. 



2.2.1 The quark current A B 
The quark current is given by 

^ Q =*7 M 75* (16) 

where \I/ should be understood to be the bagged quark field. Therefore the quark current con- 
tribution to the FSAC is given by 

o° BQ = (p\J B d 3 r^ 3 ^\p). (17) 

The calculation of this type of matrix elements in the CBM is nontrivial due to the baryon 
charge leakage between the interior and the exterior through the Dirac sea. But we know how 
to do this in an unambiguous way. A complete account of such calculations can be found in 



|3|, 16, |15|| . The leakage produces an R dependence which would otherwise not be there in the 
matrix element of Eq.(|l7|), as shown in Fig. 1. It is significant that as seen in the figure there 
is no contribution for zero radius, that is in the pure skyrmion scenario for the proton. The 
contribution grows as a function of R towards the pure MIT result that would technically be 
reached for infinite radius. The result of this calculation was first presented in refs. [||, No 
new ingredient has been added. 

2.2.2 The meson current A% 

Since we shall not add anything new to our previous result obtained in H, ||], we will just quote 
the result. Due to the coupling of the quark and rj fields at the surface, we can simply write the 
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rj contribution in terms of the quark contribution, 

< = ,nl + U 2 <Pl / d 3 r^ 737 5^|p). (18) 
' 2(1 + y v ) + y^ J B 

where y^ = m^R. In Fig. 1 we show the radial dependence of this contribution, which arises 
from the charge leakage mechanism, and follows the quark distribution. Since the rj field has 
no topological structure, its contribution also vanishes in the skyrmion limit. This illustrates 
how the dynamics of the exterior can be mapped to that of the interior by boundary conditions. 
We may summarize the analysis of these two contributions by stating that no trace of the CCP 
is apparent in Fig. 1. Thus if the CCP were to emerge, the only possibility would be that the 
gluons do the miracle! 

2.2.3 The gluon current A Bq 

Understanding the FSAC and its implications in the present framework involves crucially the 
role of the gluon contribution, in particular its static properties and vacuum fluctuations, i.e., 
the Casimir effects. The calculation of the Casimir effects constitutes the principal aim of this 
work. 

We begin by dividing the gluon current into two pieces 

Ab g = ^G,stat + ^G,vac- ( 19 ) 

The first term arises from the quark and r/ sources, while the latter is associated with the 
properties of the vacuum of the model. One might worry that this contribution could not be 
split in these two terms without double counting. That there is no cause for worry can be seen 
in several different ways. Technically, it is easy to check it by noticing that the former acts on 
the quark Fock space and the latter on the gluon vacuum. Thus, one can interprete the former 
as a one gluon exchane correction to the quanity. One can also show this intuitively by making 



the analogy to the condensate expansion in QCD [25], where the perturbative terms and the 
vacuum condensates enter additively to the lowest order. 

Let us first describe the static term. We assume initially for simplicity that there is no rj 
coupling. Then the boundary conditions for the gluon field would correspond to the original 
MIT ones [p6fl . The quark current is the source term that remains in the equations of motion 
after performing a perturbative expansion in the QCD coupling constant, i.e., the quark color 
current 

9*07/^*0 (20) 
where the fields represent the lowest cavity modes. In this lowest mode approximation, the 
color electric and magnetic fields are given by 

A a r 
' 47r r 2 



Ea =9s-- 2 p(r) (21) 



§a = (t? (3 ^ • f - *> + ^ + 2M(r))(7 ) (22) 

where p is related to the quark density p' asf] 

p(r,T) = f dsp'(s) (23) 



2 Note that the quark density that figures here is associated with the color charge, not with the quark number 
(or rather the baryon charge) that leaks due to the hedgehog pion. 
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and fj,, M to the vector current density 



fi(r) = / dsfi'(s 



o 



M{r) = J ds 



./Af) 

s 3 ' 



The lower limit V is taken to be zero in the MIT bag model - in which case the boundary 
condition is satisfied only globally, that is, after averaging - and r = R in the so called monopole 
solution ^] - in which case, the boundary condition is satisfied locally. 

After having clarified the procedure by recalling our old calculation and pointing out the 
main difference of the present with respect to that one, we proceed to introduce the rj field. 
We perform the same calculation with however the color boundary conditions Eqs.(j2[) and (||) 
taken into account. In the approximation of keeping the lowest non-trivial term, the boundary 
conditions become 

r-^ tat = -^r-B^r,(R) (24) 

r x B a stat = x E«r)(R). (25) 

Here Eg and Bg are the lowest order fields (f| given by (|2l| ) and ( |22|) and rj(R) is the meson 
field at the boundary. The rj field is given by 

where the coupling constant is determined from the surface conditions || 0. 

Note that the magnetic field is not affected by the new boundary conditions, since Eg points 
into the radial direction. The effect on the electric field is just a change in the charge, i.e., 



where 



«<*(!-) =p(r,r)+p,(B) (27) 



The contribution to the FSAC arising from these fields is determined from the expectation 
value of the anomaly 

a° G ,stat = (P\ ~ I d 3 rx 3 E^ tat ■ B a stat \p). (29) 

vr Jb 

The result of this contribution is shown in Fig. 2, where we show the MIT solution, the monopole 
one and the correction associated to both due to the color couplingP]. One sees that including 
the n contribution in p s tat(f) brings a non-negligible modification to the FSAC but does not 
modify the result qualitatively. What is most striking is the drastic difference between the 
effect of the MIT-like electric field and that of the monopole-like electric field: The former 
is totally incompatible with the Cheshire Cat property whereas the latter remains consistent 
independently of whether or not the rj contribution is included in p s tat- 

3 We have also investigated electric fields of the form {-^ + Br)f, but the results do not change much with 
respect to the ones shown since the B term tends to be small. 
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Figure 2: Dependence of a° G stat on the choice of V and the boundary conditions as a function 
of bag radius : (a) with an MITdike electric field without r] coupling, (b) with a monopole-like 
electric field without r\ coupling, (c) with an MITdike electric field with r] coupling, and (d) with 
a monopole-like electric field with r\ coupling. 
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Finally the Ac, vac term arises from the so called Casimir effect of the anomaly term. The 
vacuum in the cavity and the perturbative vacuum are different due to the geometry of the 
cavity. This effect has been considered for many other observables and also for the quarks in 
this calculation [JO], 15, Q, but never for the gluons. We proceed in the next section to describe 



this calculation. 



3 The Gluon Casimir Calculation 

The quantity that we wish to calculate is the gluon vacuum contribution to the flavor singlet 
axial current of the proton. It can be done by evaluating the expectation value 

( 0b | - ^± [ d 3 rx 3 (E a • B a )\0 B ) (30) 

where \0b) denotes the vacuum in the bag. The standard way to evaluate this expectation 
value would be to expand the field operators in terms of the classical eigenmodes that satisfy 
the equations of motion and the boundary conditions. Although well-defined, this approach is 
technically involved. We have not yet obtained any quantitative results to report. In this paper, 
we shall proceed in the opposite direction. Instead of arriving at the CCP as in the standard 
approach, we shall assume the CCP and evaluate the Casimir contribution with the expression 
that follows from the assumption. The idea goes as follows. 

The CCP states that at low energy, hadronic phenomena do not discriminate between QCD 
degrees of freedom (quarks and gluons) on the one hand and meson degrees of freedom (pions, 
etas,...) on the other, provided that all necessary quantum effects (e.g., quantum anomalies) are 
properly taken into account. If we consider the limit where the rj excitation is a long wavelength 
oscillation of zero frequency, the CCP asserts that it does not matter whether we choose to 
describe the 77, in the interior of the infinitesimal bag, in terms of quarks and gluons or in terms 
of mesonic degrees of freedom. This statement, together with the color boundary conditions, 
leads to an extremely simple and useful local formula 



E B ~— ^r- f - 2 G, (31) 

where only the term up to the first order in 77 is retained in the right-hand side. Here we adapt 
this formula to the CBM. This means that the couplings are to be understood as the average 
bag couplings and the gluon fields are to be expressed in the cavity vacuum through a mode 
expansion. In fact, by comparing the expression for the rj mass derived in using Eq.(|3l|) 
with that obtained by Novikov et al [^] in QCD sum-rule method, we note that the matrix 
element of the G 2 in ( |3l| ) should be evaluated in the absence of light quarks. This means, in 
the bag model, the cavity vacuum. That the surface boundary condition can be interpreted 
as a local operator is a rather strong CCP assumption which while justifiable for small bag 
radius, can only be validated a posteriori by the consistency of the result. This procedure is the 
substitute to the condensates in the conventional discussion. 



Substituting Eq.(31) into Eq.fl30|) we obtain 



(0 B |-^£ / al 3 rx 3 (E a ■ B a )\0 B ) 

TT JV 
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^-(p\S 3 \p)(0b\ I d\ l -G 2 x z x z \Q B ) 
Jv I 

(p\S 3 \p)(N*-l) 

£ I d 3 r(K -En-En- E n )x 3 x 3 , (32) 

where we have used that rj has a structure of (S ■ r)y(R). Since we are interested only in the 
first order perturbation, the field operator can be expanded by using MIT bag eigenmodes (the 
zeroth order solution). Thus, the summation runs over all the classical MIT bag eigenmodes. 
The factor (iV c 2 - 1) comes from the sum over the abelianized gluons. 



The next steps are the numerical calculations to evaluate the mode sum appearing in Eq.(32): 

(i) introduction of the heat kernel regularization factor to classify the divergences appearing in 
the sum and (ii) subtraction of the ultraviolet divergences. 

3.1 Normalization of the eigenmodes 

The classical eigenmode of the (abelianized) gluons confined in the MIT bag can be classified 
by the total spin quantum numbers (J, M) given by the vector sum of the orbital angular 
momentum L and the spin S, 

J = L + S, (33) 

and the radial quantum number n. There are two kinds of classical eigenmodes according to the 
relations between the parity and the total spin; (i) M-mode with the parity it = and 

(ii) E-mode with the parity ir = — (— 1) J+1 . Here, the extra minus sign is due to the negative 
intrinsic parity of gluon. 

We will work with the vector fields with the gauge choice, 

G = 0, and V • G = 0. (34) 
Then the electric field and the magnetic field are obtained through the relations 

S = -f , (35) 
B = V x G. (36) 
Explicitly, the solutions are obtained as 
(i) M-modes : 



G 



kU)W = ^MjjMY W/I (r), (37) 



(ii) E-modes 



jj+i(unr)Yjj +1<M (f) (38) 



2J + 1 



+ V 2J+ 1 l Jj ~ 1 ^ wr ^ J ' J ~ 1 ' M ^^ 



(39) 
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where Yjj m is the vector spherical harmonics of the total spin J composed of the angular 
momentum i and je(x) is the spherical Bessel functions. The energy eigenvalues are determined 
to satisfy the MIT boundary conditions as 

(i) M-modes : 

X n j'j(X n )+jj(X n )=0, (40) 

(ii) E-modes : 

jj(Xn) = (41) 

where we have defined X n = to n R. The normalization constants Mm,e will be specified below. 
The field operator G(f, t) is expanded in terms of the classical eigenmodes as 

G(r,t) = £ (a M G M (r>— 4 + a\ v} G\ u} {r)e +i ^) , (42) 
M 

where {z^} denotes the quantum number set (n, J, M, A =E or M). 

We determine the normalization constants Mm E hi such a way that the free gluon Hamilto- 
nian operator 

If, , , , , (, 3 ) 



H = - I d 3 r(E -E + B-B) 
2 



becomes 



when Eq.(42) is substituted into Eq. 
eigenmodes as 



(44) 

M 

It leads to a normalization condition for the classical 



1 



Then the normalization constants are determined explicitly as 



M e 



X n R 

X n R 2 j'j_i(X n ^ 



-1/2 



(45) 

(46) 
(47) 



3.2 Matrix elements 

The first step is to calculate the matrix elements 

73 



Q{u} = J d A r{B\ v} ■ B {v} - E\ v} ■ E {u} )x 3 x 3 . 



^From Eq.(|37j), we obtain 

E {u}(f) = (+iu n )MMjj(vnr)Yj !Jt M(r), 



B {u}{r) = (+iLU n )M M 



J 



2J + 1 



jj +1 (aj n r)Yj^ +ltM (r) 



+y 2j + + 1 3J-i( w n r )Yj,J-iM f ) 



(48) 

(49) 
(50) 

(51) 
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Figure 3: jj(x) and I(x) as a function of x. 

for the M-modes and the similar equations with E and B being interchanged for the E-modes. 
We encounter in the calculation the following angular integrals 



<MY* JAM -Yj AM xl. (52) 

By using that x| = (4/3)^7r/5l20 + 1/3 and the Wigner-Eckart theorem, we obtain 

J dQYj i M ■ Yj AM xl = cj4 J( J + 1) - 3M 2 ) + ^ (53) 

where cj^ is a constant that depends only on J and £. We have to perform the summation over 
M, which runs from — J to J, which cancels the contribution of the first term, therefore we can 
take effectively 1/3 as the result of the integral. 

Finally, we obtain the matrix elements for the M-modes as 

n (M) _ 1 x3dx ~ 27TI-7j+i( a; ) ~ ^T+T^J-i^) . . 



3 - ij-i(X n ) JJ+1 (X n )] 

In the case of the E-mode, we obtain exactly the same formula except the minus sign in front 
of it. (Note that the formulas for the electric field and the magnetic field are interchanged.) 

We have found that the matrix elements for the E-mode vanish up to our numerical accuracy 
as shown in Fig. 3. Here, the solid line is the spherical Bessel function jj{x) and the dashed line 
is the integral 

f x r j j + 1 

H x ) = J o V 3 dy jj(y) - 2J + - 2J + ^ j-iiv) 



(55) 



We see that the zeroes of I(x) and j(x) coincide, thus showing that Q n \X n ) = 0. We have 
been unable however to prove this result analytically, except the trivial case of J = 0. 

3.3 The mode sum 

In order to regularize the mode sum, we introduce a heat kernel factor exp(— rX n ); 

S(r)^J2(2J + l)Q^e-^, (56) 

n,J 
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Figure 4: Diverging properties of S(t) as a function of the heat kernel regularization parameter 
r. All the magnetic modes up to u n R(= X„)=100(solid circle), 150(solid square), 200(solid 
diamond) and 250(solid triangle) are included in the sum. 



where we have carried out the trivial sum over M and the vanishing E-mode contribution is 
excluded. 

Fig. 4 shows the numerical results of the sum up to X max =100, 150, 200, 250 for the 40 values 
of r from 0.0025 to 0.1 with the step 0.0025. We can see that below r < 0.06 the convergence is 
poor. However, it is enough to see the presence of an 1/t 2 divergence. If we fit the data above 
r > 0.06, we obtain 

s(T) = aio6i_ao8i6 + 0478 _ 00285T 

Apart from a possible logarithmic divergence, there are quadratic and linear divergences as we 
set r equal to zero . We shall remove these divergences following a procedure commonly used in 



Casimir problems [29]. Caveat on this procedure will be highlighted in the discussion section. 
Now if we neglect logarithmic divergences that might be present , the best way to get rid of the 
quadratic and linear divergences is to evaluate 

S(t) + 2tS'(t) + \t 2 S"{t) = ^ (2 J + l)Q„,j(l - 2rX n + 0.5T 2 X 2 ) e - rX " . (58) 

~~ n,J 

Fig. 5 show the results on this quantity for 80 values of r ranging from 0.0025 to 1. We 
see that no serious divergences appear anymore. By fitting the convergent data with the above 
expressions for r, we obtain for the finite part of the sum 0.0478, from the cubic function fit, 
and 0.0456, from the quadratic one. These results are comparable to the finite term of the above 



naive fitting procedure (57), which yielded 0.0478 



Once we have the numerical value on the mode sum, the gluon vacuum contribution to FSAC 
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Figure 5: 5(r) — 2tS'(t) + ^r 2 5"'(r) as a function of r. The finite term of 5(r) is extracted by 
fitting these quantities to a cubic and quadratic curves. 



can be evaluated simply as 



u G,vac 



where y{R) is related to 

y(R) = - 



(2.10) 2 8 y (R) , , 

^ '— x - x yv 7 x 0.0478 , 

2 2 122MeV v h 



(59) 



3m 2 



(1 + m v R) 



brfn [2(1 + m^i?) + (m v R) 2 ][m n R) 



2 U B Q - 



(60) 



We have used N F = N c = 3, a s = 2.2, /„ = y/N^/2f^ ~ 122MeV and m,, = 958 MeV. 



4 Results and Conclusions 

Our numerical results are given in Fig.^. Standard MIT bag parameters were used for the 
calculation. The quarkish component of the FSAC is given by the sum of the quark and rj 
contributions, ag Q +a® and the gluonic component by a° G stat +a°c vac - Both increase individually 
as the confinement size R is increased but the sum remains small, < a^ ota i < 0.3 for the whole 
range of radii, consistent with the experiment, a«* = a°(oo) = 0.10i°;{£. It is remarkable that 
a{R = 0) ~ a{R ~ 1.5 fm) while each component can differ widely for the two extreme radii. 

We have shown that the principal agent for the observed small FSAC in the proton in the 
framework of the chiral bag model is the Cheshire Cat phenomenon (CCP). It is the CCP that 
assures the cancellation between two contributions, one from the quarkish component and the 
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Figure 6: Various contributions to the flavor singlet axial current of the proton as a function 
of bag radius and comparison with the experiment: (a) quark plus r\ contribution (a% + a^), 
(b) the contribution of the static gluons due to quark source 

stat)' ( c ) gluon vacuum 

contribution (aQ vac ), and (d) their sum (a® otal ). The shaded area corresponds to the range 
admitted by experiments. 



other from the gluonic component in the particular way the separation is made. For a small bag 
radius, both components are small, so the net FSAC being small is immediate. This is consistent 
with the observation that in the limit that R — > 0, we recover the skyrmion description which 
gives a vanishing FSAC at the leading order, modified by matter fields at the next order. At 
large bag radius which leads to the MIT bag model, both the quarkish contribution and the 
gluonic contribution are large but they cancel. Our assertion is that this cancellation is caused 
by the CCP. We should however recall that the separation between the quarkish component 
and the gluonic component we adopted in (14) and ( |T5| ) is entirely arbitrary although the sum 
is unique. Whether the separate component by itself is large or small has no physical meaning. 
Only the total does. Different separations would lead to different scenarios leading to the same 
small value. It is plausible that in some limit - unknown to us - the FSAC would be exactly 
zero with the finite nonzero value indicating a departure from this limit. Understanding this 
limit would allow a unique separation of the components. 

One of the principal results of this paper is that it is possible to have a nonzero value for the 
FSAC at R = and is of the same size as at large R^. While the effect of the surface color 
anomaly term is generally small for all radii, the finite nonzero value of FSAC for R = is 
assured by the surface boundary term. Thus the violation of the CCP observed in the previous 
calculations at R = || 0] is neatly eliminated by the color anomaly boundary condition. More 
importantly, the monopole structure of the color electric field previously proposed is found to be 
required for the sign that comes with the important static gluonic contribution from the quark 



4 The reason for this nonzero value is intimately connected with the CCP, since it is the finite part of the gluon 
mode sum which normalizes the value of this contribution at the origin. Moreover the color boundary condition 
provides us with a decreasing r\ field contribution which changes softly as a function of R. 
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source. We believe that this cancellation is a manifestation in the bag scenario of the recently 
discovered one for QCD [pi]]. The MIT configuration would strongly violate the CCP. We are 
thus led to the conclusion that the CCP requires the monopole configuration for the color electric 
field. Whether or not this configuration leaves undisturbed other - successful - phenomenology 
was discussed in || . 

In calculating the gluonic Casimir effect, we made the ah initio assumption that the CCP 
holds, an assumption which is expected to be valid for small bag radius. We then extend it, 
in accordance with the CCP, to all bag radii. We can justify this only a posteriori by showing 
that the CCP assumption is consistent with what one gets out. Note however that the gluonic 
Casimir effect is most significant for small R where it is needed for the CCP and plays little 
role for large R. Thus our assumption is validated. It would of course be more satisfying if one 
could obtain the CCP as an output of the formalism, not put in as an input. Such a calculation 
is in progress. 

We should mention a caveat left unspecified in the text in regularizing this Casimir contri- 
bution. Since a G vac vanishes when the rj field is removed, the so-called "vacuum contribution" 
is duly subtracted in what we have computed. However we have also explicitly subtracted 
quadratic and linear divergences appearing from the mode sum by resorting to a procedure used 
in the past in most of Casimir- type calculations [29] which as far as we know, is physically rea- 
sonable but has not yet been rigorously justified from first principles. The same caveat applies 
to our calculation as it does to others. The finite term we have obtained might therefore be 
subject to additional finite corrections. In this paper we have invoked the Cheshire Cat Principle 
to ignore such corrections in a Gvac . We hope that the calculation in progress [18] will eliminate 
this ambiguity. 

Given the caveat mentioned above and the approximations used, our result can at best be 
qualitative. A better treatment (such as a more realistic gauge coupling constant running with 
the bag size, a more accurate calculation of a Gvac etc.) might modify the result quantitatively. 

Even so, we believe it to be quite robust that the overall FSAC is small, ~ 0.3 and that it is 
more or less independent of the confinement size. 



Acknowledgments 

This work started as a pilot activity in astro-hadron physics at Korea Institute for Advanced 
Study (KIAS) when two of us (MR & VV) were visiting the Institute in February 1998. They 
would like to thank KIAS for support and C W. Kim, Director of KIAS, for his efforts to make 
the collaboration a most profitable and pleasant one. VV also benefited from a travel grant 
from the Conselleria de Cultura, Educacio i Ciencia de la Generalitat Valenciana, which allowed 
his visit to DAPNIA-CEA Saclay, where part of this work was done. He thanks J. P. Guichon 
and J.M. Laget for their efforts to make the visit possible and the members of the Service for 
the hospitality. His work was partially supported by DGICYT-PB97-1227 and by the TMR 
programme of the European Commission ERB FMRX-CT96-008. The work of HJL, DPM and 
BYP is partially supported by the Korea Science and Engineering Foundation through CTP of 
Seoul National University. 



16 



References 

[1] H.B. Nielsen and A. Wirzba, Elementary Structure of Matter, Springer Proceedings in 
Physics 26 (1988); Proc. on the fundamental structure of hadronic matter, Les Houches 
France (1987); S. Nadkarni, H.B. Nielsen and I. Zahed, Nucl. Phys. B253 (1985) 308; S. 
Nadkarni and H.B. Nielsen, Nucl. Phys. B263 (1986) 1; S.Nadkarni and I. Zahed, Nucl. 
Phys. B263 (1986) 23; R. Perry and M. Rho, Phys. Rev. D34 (1986) 1169. 

[2] For detailed review, see, e.g., Chiral Nuclear Dynamics by M.A. Nowak, M. Rho and I. 
Zahed (World Scientific Pub. Singapore, 1996) Chapter 8. 

[3] S. Coleman, Phys. Rev. D 11 (1975) 2088; S. Mandelstam, Phys. Rev. D 11 (1975) 3026; 
A. Ferrando and V. Vento, Phys. Rev. D 49 (1994) 3044. 

[4] S. Weinberg, Physica 96A (1979) 327; The Quantum Theory of Fields //(Cambridge Press, 
1996); Nature 386 (1997) 234; H. Georgi, Weak Interactions and Modern Particle Physics 
(The Benjamin/Cummings Pub. Co., Menlo Park, CA, 1984). 

[5] B.-Y. Park, V. Vento, M. Rho and G.E. Brown, Nucl. Phys. A504 (1989) 829. 

[6] B.-Y. Park and V. Vento, Nucl. Phys. A513 (1990) 413. 

[7] M. Rho and V. Vento, Nucl. Phys. A 662 (1997) 538. 

[8] F. Yndurain, The theory of quark and gluon interactions 2nd ed . (Springer Verlag, Heidel- 
berg 1993); T.Mutha, Foundations of Quantum Chromodynamics (World Scientific, Singa- 
pore 1987). 

[9] J. Gasser and H. Leutwyler, Ann. Phys. 158 (1984) 142; T.H.R. Skyrme, Nucl. Phys. 31 
(1962) 556; E. Witten, Nucl. Phys. B223 (1983) 422; ibid. 433. 

[10] V. Vento, M. Rho, E.M. Nyman, J.H. Jun and G.E. Brown, Nucl. Phys. A345 (1980) 413; 
G.E. Brown, A.D. Jackson, M. Rho and V. Vento, Phys. Lett. B140 (1984) 285. 

[11] R.D. Ball, Perturbative Evolution of Polarized Structure Functions, Erice 1995, Eds. B. 
Frois, VW. Hughes and N. de Groot (World Scientific, Singapore 1997). 

[12] J.Ellis and R.L.Jaffe, Phys. Rev. D9 (1974) 1444. 

[13] G. 't Hooft, Phys. Rep. 142 (1986) 357. 

[14] M. Rho, A.S. Goldhaber, G.E. Brown, Phys. Rev. Lett. 51 (1983) 747; J. Goldstone and 
R.L. Jaffe, Phys. Rev. Lett. 51 (1983) 1518; G.E. Brown, A.D. Jackson, M. Rho and V. 
Vento, Phys. Lett. B140 (1984) 285; A.D. Jackson, L. Vepstas, E. Wrist and D. Kahana, 
Nucl. Phys. A462 (1987) 661. 

[15] B.-Y. Park and M. Rho, Z. Phys. A331 (1988) 151. 

[16] B.-Y. Park, Z. Phys. A331 (1988) 473. 

[17] H.B. Nielsen, M. Rho, A. Wirzba and I. Zahed, Phys. Lett. B269 (1991) 389. 



17 



[18 
[19 

[20 

[21 
[22 

[23 
[24 

[25 

[26 
[27; 
[28 



H.-J. Lee, D.-P. Min, B.-Y. Park, M. Rho and V. Vento, work in progress. 

S. Forte, Proceedings of PANIC96, Williamsburg, May 1996 and SPIN96, Amsterdam, 
September 1996, |hep-ph/960l"238 . 

For summary, J. Ellis and M. Karliner, "The Strange Spin of the Nucleon," Lectures at the 
Int. School of Nucleon Spin Structure, Erice, August 1995, hep-ph/960128C| , 



N. Kochelev and V. Vento, jhep-ph/9809238 



G.M. Shore and G. Veneziano, Phys. Lett. B244 (1990) 75, Nucl. Phys. B381 (1992) 23; 
|hep-ph/97092T3 , 



G. Altarelli, R.D. Ball, S. Forte and G. Ridolfi, Acta Phys.Polon. B29 (1998) 1145. 

S.B Treiman, R. W. Jackiw, B. Zumino, and E. Witten in Current algebra and anomalies 
(World Scientific, Singapore, 1996). 

M.A. Shifman in Vacuum structure and QCD sum rules (North Holland, Amsterdam 1992); 
S. Narison in QCD spectral sum rules (World Scientific, Singapore 1989). 

T. De Grand, R.L. Jaffe, K. Johnson and J. Kiskis, Phys. Rev. D12 (1975) 2060. 

H. B. Nielsen, M. Rho, A. Wirzba and I. Zahed, Phys. Lett. B281 (1992) 345. 

V. A. Novikov, M.A. Shifman, A.I. Vainstein and V.I. Zakharov, Nucl. Phys. B 191 (1981) 
301. 



[29] J. Schwinger, Lett. Math. Phys. 24 (1992) 59, 227; Proc. Natl. Acad. Sci USA 89 (1992 ) 
4091, 11118; K.A. Milton and Y. Jack Ng, Phys. Rev. E57 (1998) 5504, [hep-th/97071~22l 



18 



